#p = ms.readMS( 'msOuts/out10_10000.ms' )
#
#mySFS = SiteFrequencySpectrums(p)
#eSFS = mySFS.getExpectedSFS()
#
#print eSFS
#print
#
#twoSidedSFS = mySFS.getTwoSidedSFS( 0.10 )
#print twoSidedSFS
#print
#
#myESFS = np.ravel(eSFS[0, 1:])
#lower = np.ravel( twoSidedSFS[0, 1:] )
#upper = np.ravel( twoSidedSFS[1, 1:] )
#
#sfsX = range(1, myESFS.size+1)
#
#plt.plot( sfsX, myESFS, 'b', sfsX, lower, 'r--', sfsX, upper, 'r--' )
#plt.axis( [.9, myESFS.size + 0.1, -0.1, 1.1] )
#plt.ylabel('Frequency')
#plt.xlabel('n-tons')
#plt.title('Expected SFS with 0.1 and 0.9 levels, n=10')
#plt.show()
#        
